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Abstract 

Measuring the T(h — > 77)Br(/i — * bb) decay at the photon collider at TESLA is studied 
for a Standard Model Higgs boson of mass = 120 GeV. The main background due 
to the process 77 — ► QQ(g), where Q = b, c, is estimated using the NLO QCD program 
(G. Jikia); the results obtained are compared with the corresponding LO estimate. Using 
a realistic luminosity spectrum and performing a detector simulation with the SIMDET 
program, we find that the T(h — > 77)Br(/i — > bb) decay can be measured with an accuracy 
better than 2% after one year of photon collider running. 
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1 Introduction 



A search of the last missing member of the Standard Model (SM) family, the Higgs boson, is 
among the most important tasks for the present and future colliders. Once the Higgs boson is 
discovered, it will be crucial to determine its properties with a high accuracy. A photon collider 
option of the TESLA e + e~ collider [1] offers a unique possibility to produce the Higgs boson as 
an s-channel resonance. The neutral Higgs boson couples to the photons through a loop with 
the massive charged particles. This loop-induced /177 coupling is sensitive to contributions of 
new particles, which appear in various extensions of the SM. 

The SM Higgs boson with a mass below ~ 140 GeV is expected to decay predominantly 
into the bb final state. Here we consider the process 77 — > h — > bb for a Higgs-boson mass of 
m h = 120 GeV at a photon collider at TESLA. Both the signal and background events are 
generated according to a realistic photon-photon luminosity spectrum [2], parametrized by a 
CompAZ model [3]. Our analysis incorporates a simulation of the detector response according 
to the program SIMDET [4]. 



2 Photon— photon luminosity spectrum 

The Compton backscattering of a laser light off high-energy electron beams is considered as 
a source of high energy, highly polarized photon beams [5]. A simulation of the realistic 77 
luminosity spectra for the photon collider at TESLA, taking into account non-linear corrections 
and higher order QED processes, has become available recently [2]. In this simulation, according 
to the current design [1], the energy of the laser photons u>l is assumed to be fixed for all 
considered electron-beam energies. 

In the analysis we use the CompAZ parametrization [3] of the spectrum [2] to generate 
energies of the colliding photons. We assume that the energy of primary electrons can be 
adjusted in order to enhance the signal. Our signal production of a scalar particle corresponds 
to the case where projection of the total 77 angular momentum on a collision (z) axis J z is 
equal to zero. For y/s^ = 2E e = 210 GeV, we obtain a peak of the J z = component of the 
photon-photon luminosity spectra at the invariant mass of the two colliding photons W 77 equal 
to the considered mass of the Higgs boson, i.e. m h = 120 GeV. 

The luminosity spectra are shown in Fig. 1. The lowest invariant mass of the two colliding 
photons used in our generation of events is equal to W lini i n = 80 GeV. For the assumed ^/s^ 
value, the maximum invariant mass for the colliding photons, each produced in a single Compton 
scattering, is equal to W max i = 131.2 GeV. However, there is also a small contribution from 
the events, which correspond to the interaction of an initial electron with two laser photons 
(higher order effect). This gives a higher maximal invariant mass of the produced energetic 
photon beams, namely H / max 2 = 161.5 GeV. The results presented in this paper were obtained 
for an integrated luminosity of the primary e~e~ beams equal to L^ om = 502 fb^ 1 , expected 
for one year of the photon collider running [2] . The resulting 77 luminosity is then expected to 
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Figure 1: Photon-photon luminosity spectra for y^s^T = 210 GeV, obtained with CompAZ 
parametrization of Telnov's simulation, as a function of the invariant mass of two colliding 
photons W 77 . The contributions of states with the total 77 angular momentum projected on a 
collision (z) axis, J z =0 and J z = ±2 (denoted as J z — 2), are shown separately. 

be: L 77 = 409 fb _1 , or 84 fb" 1 for W, n > 80 GeV . 

In the earlier analyses, for instance in [6], the spectra that were used had been derived 
from the lowest order QED calculation for the Compton scattering, with a fixed parameter 
x — 4E m2 L equal to 4.8. The realistic spectrum [2], parametrized by the CompAZ model, differs 
significantly from the spectrum of the high-energy photons used in [6], which is shown in Fig. 2. 
The comparison is made for two combinations of the helicities of two colliding photons, (±, ±) 
and (±, =f), with the total angular momentum projected on a collision (z) axis equal to and 
±2, respectively. 



3 Details of a simulation and the first results 

We calculate the total width and branching ratios of the SM Higgs boson, using the program 
HDECAY [7], where higher order QCD corrections are included. A generation of events was 
done with the PYTHIA 6.205 program [8], with the parameters for a Higgs boson as in the 
HDECAY. A parton shower algorithm, implemented in PYTHIA, was used to generate the 
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Figure 2: Photon-photon luminosity spectra used in the analysis of the SM Higgs-boson pro- 
duction with mass = 120 GeV, as a function of the invariant mass of two colliding photons 
W 77 . The spectrum used here, as obtained from CompAZ parametrization based on Telnov 
simulation (hatched areas), is compared with a spectrum derived from the lowest order QED 
predictions for the Compton scattering, used in the earlier analysis (lines). The total luminosity 
distribution ( J z = 0, ±2) and the J z = contribution are shown, separately. 



final-state particles. 

The background events due to processes 77 — > bb(g), cc(g) were generated using the pro- 
gram written by G. Jikia [6], where a complete NLO QCD calculation for the production of 
massive quarks is performed within the massive-quark scheme. The program includes exact 
one-loop QCD corrections to the lowest order (LO) process 77 — > bb, cc [9], and in addition 
the non-Sudakov form factor in the double-logarithmic approximation, calculated up to four 
loops [10]. 

For a comparison we generate also the LO background events, using the QED Born cross 
section for the processes 77 — > bb and 77 — > cc, including in addition a parton shower, as 
implemented in PYTHIA. 1 

The fragmentation into hadrons was performed using the PYTHIA program. A fast simula- 
tion for a TESLA detector, the program SIMDET version 3.01 [4], was used to model a detector 

-"^For consistency with Jikia's program, we use a fixed electromagnetic coupling constant equal to a em 
1/137. 
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performance. The jets were reconstructed using the Durham algorithm, with y cut = 0.02; the 
distance measure was defined as yij = 2 min(Ef, Ej){\ — cos^)/ E% is , where E vis is defined as 
the total energy measured in the detector. 

The double fo-tag was required to select the signal h — > bb events. Since no suitable flavour- 
tagging package exists for the SIMDET 3.01 program 2 , we assume, following the approach used 
in [6], a fixed efficiency for the 66-tagging, equal to e bb = 70%, and a fixed probability for a 
mistagging of the cc events, a main background to the bb events, equal to e cc = 3.5%. 

The following cuts were used to select properly reconstructed bb events: 

1. a total visible energy E vis greater than 90 GeV, 

2. since the Higgs boson is expected to be produced almost at rest, the ratio of the total 
longitudinal momentum of all observed particles to the total visible energy is taken to be 
\P Z \/E ms < 0.1, 

3. a number of jets Nj ets = 2, 3, so that events with one additional jet due to hard-gluon 
emission are accepted, 

4. for each jet, i — 1, ...,Nj ets , we require | cos^| < 0.75, where cos#j = p Z i/\pt\- 

Using the above cuts, we obtain the distributions of the reconstructed 77 invariant mass W rec 
for a signal and for a background, shown in Fig. 3. The J z = and J z = ±2 contributions 
from the NLO background, with bb(g) and with cc(g) final states, are shown separately. For a 
comparison, the estimated LO background is presented as well (dotted line). Note that the NLO 
background contribution is approximately two times larger than the LO one. This is mainly 
due to the J z = component, which is strongly suppressed in the LO case; however in the case 
of NLO it gives a large contribution, especially pronounced in the high-W 7 ^ part of the signal 
peak. A more detailed comparison of the LO and NLO background estimations is presented in 
Appendix. Other background contributions, from the resolved photon(s) interactions and the 
overlaying events, were found to be negligible [13]. 



4 Final results 

Assuming that the signal for Higgs-boson production will be extracted by counting the bb events 
in the mass window around the peak, and subtracting the background events expected in this 
window, we can calculate the expected relative statistical error for the partial width multiplied 
by the branching ratio, T(h — * r yy)Bi(h — > bb), in the following way: 



A 



F(h -> 7 7 )Br(> -> bb) 



obs 



T(h -> 77)Br(/i -> bb) 



N obs - Ni 



bkgd 



2 The 6-tagging code adapted for the new version of SIMDET 4.01 [11] should become available soon [12]. 
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Figure 3: Reconstructed invariant mass W rec distributions for the selected bb events. Contribu- 
tions of the signal, due to the Higgs boson with a mass = 120 GeV, and of the heavy-quark 
background, calculated in the NLO QCD, are indicated. For comparison, the LO background 
estimate is also plotted (dots). Arrows indicate the mass window optimized for the measure- 
ment of the T(h — > 77) Br (ft, — ► bb). 

The accuracy expected for the considered quantity T(h —>■ 77)Br(ft — ► bb), if estimated from 
the reconstructed invariant-mass distribution obtained for the Higgs-boson mass of 120 GeV 
in the selected mass region, between 106 and 126 GeV (see Fig. 3), is equal to 1.9%. It is in 
agreement with the result of a previous analysis [6]. 

A long tail in the reconstructed mass W rec distribution obtained for the h — ► bb events seen 
in Fig. 3 is due to the escaping neutrinos, which mainly originate in the semileptonic decays of 
the D- and I?-mesons. This tail can be effectively suppressed by applying an additional cut on 
Pt/Et, where Pt and Et are the absolute values of the total transverse momentum of an event 
Pt and the total transverse energy, respectively. 3 The cut relies on demanding the Pt/Et to 
be small. The W rec distributions for the h — > bb events obtained by applying various Pt / Et 
cuts are shown in Fig. 4. In this figure we use different colours to denote the different total 
energies of neutrinos in the event, E vs . The effects due to the detector resolution influence a 
shape of the distribution for W rec > rrih, whereas for a lower W rec one observes a significant 
effect of the escaping neutrinos in the distribution. By applying a realistic Pt/Et cut, e.g. 

3 Pt (Et) is calculated as a vector (scalar) sum of the transverse momenta py. = (p%,Py,0) (the transverse 
energies E % T = E % sin#i) over all particles that belong to an event. 
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Figure 4: Reconstructed invariant mass W rec distributions for 77 — > h — > bb events, for various 
Pt/Et cuts. Contributions of events with a different total energy of neutrinos in the event E vs 
are indicated by different colours. The parameters /x and a are obtained from the Gaussian fit 
in the region (/i — a, // + 2a). 

Pt/Et < 0.04, we can obtain a mass resolution, derived from the Gaussian fit in the region 
from fi — a to + 2a, better than 2 GeV. 

Shown in Fig. 5 is the W rec distribution obtained by applying the cut P T /E T < 0.04. 
The relative accuracy expected for the T(h —>■ 77)Br(/i — > bb) measurement, calculated in the 
W rec mass range between 114 and 124 GeV (as indicated by arrows in the figure), is equal to 
2.2%. We conclude that the P T /E T cut improves a mass resolution, but worsens the statistical 
significance of the measurement. 

We have found a method, which allows an increase of a signal-to-background ratio without 
reducing the event statistics. We assume that the measured missing transverse momentum is 
due to a single neutrino emitted perpendicularly to the beam line. 4 Then, we introduce the 



4 Due to a large spread of the photon beam energy, no constraints can be imposed on the longitudinal 
momentum. 
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Figure 5: As in Fig. 3, for the reconstructed invariant mass W rec distributions for the selected 
bb events, obtained by applying an additional cut on the ratio of the total transverse momentum 
and total transverse energy, Pt/Et < 0.04. 



corrected, reconstructed invariant mass as: 



W corr = JW r \ c + 2P T (E vis + P 7 



(1) 



The distributions of the W corr , obtained for the signal and background events, are shown in 
Fig. 6. The most precise measurement of the Higgs-boson production cross section is obtained 
using the mass window W corr between 115 and 128 GeV, as indicated by arrows. In the selected 
W corr region one expects, after one year of photon collider running at nominal luminosity, 
about 5900 reconstructed signal events and 4600 background events (i.e. S/B ps 1.3). This 
corresponds to the expected relative statistical precision of the measurement: 



T(h -> 77)Br(/i -> bb) 



T(h -> 77 )Br(/i -> bb) 



1.7% 



5 Conclusions 



Our analysis shows that, for the SM Higgs boson with a mass around 120 GeV, the two-photon 
width for the bb final state can be measured in the photon collider at TESLA with a precision 
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Figure 6: As in Fig. 3, for the corrected invariant mass W corr distributions. 

better than 2%. If the reconstructed invariant mass of the event is corrected for the energy of 
escaping neutrinos from D- and 5-meson decays, we achieve a precision of the measurement 
of the T(h —>■ 77)Br(/i — > bo) equal to 1.7%. The obtained accuracy is in an agreement with 
the result of a previous analysis, based on the idealistic Compton spectrum [6] . Note, however, 
that the realistic photon-photon luminosity spectrum that we use is more challenging. 

The measurement discussed in this paper can be used to derive the partial width T(h — > 77), 
taking Br(/i — > bb) value from precise measurement at the e + e~ Linear Collider. With 1.7% 
accuracy on T(h —> 77)Br(/i — > bb), obtained in this analysis, assuming Br(/i — > bb) will be 
measured to 1.5% [14], Higgs-boson partial width T(h — > 77) can be extracted with accuracy 
of 2.3%. Using in addition the result from the e + e~ Linear Collider for Br(/t — > 77) [15], one 
can also extract r tot with precision of 10%. 

The SM Higgs-boson production 77 — > h — > bb can be considered for masses up to about 
160 GeV [6]. For higher masses of the SM Higgs boson one should consider other decay channels, 
see e.g. [16]. 
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Appendix: Comparison of the LO 

and NLO background estimates 

As it is well known, see e.g. [9, 10], the NLO corrections for the background process 77 — > bb 
are large. Whereas the LO background is strongly suppressed for a J z = contribution, this 
suppression is removed for the higher order process with an additional gluon in the final state. 
As an example, we show in Fig. 7 the ratios of the corresponding NLO and LO results for the 
W rec distribution, for the process 77 — ■> bb(g), for different values of Nj ets and J z . In each case 
one observes large differences between the NLO and LO results, both in the shape and in the 
normalization. As the precise determination of the background shape is crucial for a reliable 
estimation of the Higgs-boson width, we conclude that a rescaling of the LO estimates cannot 
be recommended as a substitute of the full NLO background analysis. 
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Figure 7: Ratio of the NLO to LO results for the reconstructed invariant mass W rec distribution 
for the background process 77 — > bb(g). 
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